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For standard interactions of neutrinos with matter bimagic baseline of length about 2540 Km 
is known to be suitable for getting good discovery limits of neutrino mass hierarchy, sin^ S13 and 
CP violation in the — > i^n oscillation channel. We discuss how even in presence of non-standard 
interactions (NSIs) of neutrinos with matter this baseline is found to be suitable for getting these 
discovery limits. This is because even in presence of NSIs one could get the i/^ ~^ ffi oscillation 
probability to be almost independent of CP violating phase 5 and 6*13 for one hierarchy and highly 
dependent on these two for the other hierarchy over certain parts of neutrino energy range. For 
another certain part of the energy range the reverse of this happens with respect to the hierarchies. 
We present the discovery limits of NSIs also in the same neutrino energy range. However, as with the 
increase of neutrino energy the NSI effect in the above oscillation probability gets relatively more 
pronounced in comparison to the vacuum oscillation parameters, so we consider higher neutrino 
energy range also for getting better discovery limits of NSIs. Analysis presented here for 2540 Km 
could also be implemented for longer bimagic baseline > 6000 Km. 
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I. INTRODUCTION 



The present experiments on neutrino oscillations confirms that there is mixing between different fiavours of 
neutrinos {u^, f^, Vt)- The probability of neutrino oscillations depends on various parameters of the neutrino 
mixing matrix-the PMNS matrix [1]. The current experiments tells us about two of the angles 6123 and 9i2 
[2] with some accuracy but for 6*13 only the upper bound is given [2] and the CP violating phase 6 is totally 
unknown. Although the mass squared difference of the different neutrinos (Am^^- = rnj — to|) are known 
to us but the sign of Aml^ (which is related to mass hierarchy) is still unknown. Due to the correlations 
among these unknowns there are ambiguities [3] in analysing neutrino oscillation datas. To reduce these 
ambiguities one may consider nciitrino oscillation experiments in long baseline [4] - particulaly in magic 
baseline [5] . The magic baseline satisfies certain condition on its length from the detector and is found to be 
about 7500 Km where the perturbative expression of probability P{i'e — i^n) becomes independent of 5 upto 
order (where a = Amli/Aml]^). Although this could result in finding out the other unknown oscillation 
parameters conveniently but for measurement of 6 this baseline is not suitable. To circumvent this problem, 
conditions on neutrino energy has been considered [6, 7] in 1^,. — > z/^ channel for which also the perturbative 
expression of probability becomes independent of S but only on a part of the neutrino energy spectrum. But 
the other part of the spectrum will be sensitive to CP violating phase 6. As pointed out in [7], one may 
consider satisfying two different energy conditions simultaneously for two hierarchies in the same baseline 
which results in fixing the length of bimagic baseline to about 2540 Km. Unlike magic baseline, here the 
baseline is shorter so the neutrino flux for such baseline is reduced by lesser amount at the detector. Also in 
this oscillation channel which has been considered in this work, the detection of muon is easier in 

comparison to some other channels where the detection of electron is required. 



We have studied the effect of non-standard interactions (NSIs) of neutrinos with matter in bimagic baseline. 
At first there is discussion on how to obtain the perturbative expression of the probability of oscillation upto 
order in Vf. — >■ channel in presence of NSIs. The NSIs present in the i/g i^i^ oscillation probability are 
Ceej ^en ^nd Eer among which Cee, SiTC ^ a but Cge has no such constraints in considering perturbation. In 
our numerical analysis we have considered the experimentally allowed range which covered the pcrtiirbative 
regime and also has gone beyond that. We have also presented the 6 and 813 independent perturbative 
expression of the oscillation probability in presence of the NSIs under two different magic energy conditions 
corresponding to two different hierarchies. In presence of NSIs different discovery limits for hierarchy of 
neutrino masses, for sin^ ^13 and also for CP violation have been shown in figures. Discovery limits of 
NSIs particularly tee, f-en and Cer for specific values of 61^ and 5 in their allowed range have also been 
presented. One may note that to satisfy the bimagic conditions one requires lower neutrino energy within 5 
GeV. However, the perturbative expression of oscillation probability shows that the NSI effect will relatively 
increase in comparison to other neutrino oscillation parameters in vacuum if the neutrino energy is higher. 
For this reason we have considered higher neutrino energy of 50 GeV also to study the discovery limits of 
NSIs in the same baseline of 2540 Km for which better limits are obtained. 
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II. z/e Vu. OCILLATION PROBABILITY IN PRESENCE OF NSI 



The fermion-neutrino interaction in matter is defined by the Lagrangian: 



Gf_ ffP\ 
V2 V'-/^' 



(1) 



where P € L,R, L = 1 — 7''', R = 1+7"'', / = e,u,d and e^^ is the deviation from standard interactions. 
There are model dependent bounds on these NSI parameters [8, 9]. In i?-parity violating Supersymmetric 
models these NSI parameters could be related to trilinear lepton number violating couplings [10]. Also such 
parameters could be sizable [11] in unified supersymmetric models [12]. The model independent bounds have 
been discussed in [13]. The above NSI parameters can be reduced to the effective parameters as: 



(2) 



where nf is the number density of fcrmion, n,. is the electron number density. In neutrino oscillation 
experiments this effective parameter (ea^) corresponds to the replacement in the matter interaction part of 
the evolution of flavoured neutrinos. This change can be seen as below: 





(3) 



In general, Cg^, Cer and e^,- 



Cer ana e^r could be complex. However, for our numerical analysis we have considered Ce^u 
and to be real. If we assumne uncorrelated errors, the bounds on can be approximately written as 
[13] 



< 



-,1/2 



L P 



(4) 



The NSIs eee, Ce/^ and eer play significant role in Ue — > Vf_i oscillation channel which we have considered. The 
bounds for these [13] in the context of neutrino oscillation for neutrinos passing through neutral earth like 
matter is Cee < 4, < 0.33 and Cei- < 3. 

In vacuum, flavor eigenstates Va may be related to mass eigenstates of neutrinos as 



Wa > -- 



Uai\i'i>, U = R23Ri3{S)Ri2 and 



1,2,3, 



(5) 



where U is PMNS matrix [1] and Rij are the rotation matrices and Ri3{5) contains the CP violating phase 
S signifying the complex rotation [11]. General probability expression for oscillation of neutrino of flavor I 
to neutrino flavor m in matter (satisfying adiabatic condition for the density of matter) is given by 



Pii^i ^ t^m) = Sim - ^Yl i?e[4f] sin^ a;, + 2 Y] Im[f^] sin 2/^^^ 



i>j 



i>j 



(6) 



where 



Jim _ TJ TJ *TT * TT 

"^ij ~ li'-' Ij '-^ mi'-' raj 1 



(7) 
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Here 

A TOy =m, - nij (9) 

and label ( ' ) indicates the neutrino matter interaction induced quantities corresponding to those quantities 
in vacuum. 

We discuss in brief the perturbative approach for evaluating the induced quantities and for obtaining the 
probability of oscillation — )• for neutrinos passing through earth matter. The diagonal neutrino mass 
matrix is approximately given by 

m w Am^idiag{0, a, 1). (10) 

The effective Hamiltonian induced by interaction of matter with neutrinos is written in weak interaction 
basis as 



H « R^sMRIs (11) 



where 



M=^Rismi^^RURMS)' + ^diag{A,0,0) + ^Rl,\x* B C \ R,s . (12) 



Z X 
X* B C 
Y* C* D, 



In equation (12) 



2Ey/2G 



A — — An? ' ^ ~ ^ ~ ^^erj Z — Aeee, B — Acfj^n, C — Ae^T, D — Ae^^j (13) 

where A is considered due to Standard model interaction of neutrinos with electron and Gf is the Fermi 
constant and rig is the electron number density written as rie = {Q.^Na)p and Na is Avogadro's number and 
p is the matter density in gm/cc. Cee, Ce/n i Cer; ^nn-, ^fir and Crr are considered due to NSI of neutrinos with 
matter. We consider magnitude of these NSI parameters except Cee not higher than a in using perturbation 
method. As Cee has been considered in the leading part of the Hamiltonian, our perturbative result is fine 
even for its highest experimentally allowed value. In our numerical analysis we have considered even the 
uppermost allowed values of other NSI parameters [13]. In equation (12), ( * ) is denoted for complex 
conjugation. 

The mixing matrix U' can be found out as U' = R23 W. Here W is the normalized eigenvectors of 
Am3iM/(2£') calculated through perturbative technique similar to the one adopted in [14]. We have taken 
into account only the non-degenerate perturbative approach. Let us consider the case where NSIs are 
present and where sin(?i3 is small and of the order of a or less. M in equation (12) can be written as 
M = M(°) + M(i) + M(2) where contains terms of the order of a\ Then we can write 

A 2 A 2 / ^^12 ^ 

MW = ^ d^agiA', 0, 1), M^^) = ^ ( b* acl, + c,,c 

\ a* 

/ s?3 -e-*^aci3s22Si3 

M(2) = ^ -e-^aci2.si2=si3 I (14) 



2E 




,-e"'aci3sf2Si3 -e*''aci2Si2Si3 -5^3 
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where 

A' = A{1 + Cee), a = C23Y + e"*''si3 + XS23, b = C23X + Ci2Q!Si2 - YS23, 

C = BC23 - C*S23, d = CC23 - DS23, e = C*C23 + BS23, f = DC23 + ^^23 (15) 

The eigenvalues of H upto second order in a are 



/ 2 
nil 

'2E 



2E 



2E 



2E 

2E 
2E 



A' + asi2 + s?3 + ^ 



\bl 

A' ' i-l + A') 



aci2 - + CC23 - idcl3 + csl^^Y 



,*|2 



'13 



{I -A') 



77 + /C23 + eS23 + (eC23 - /S23)' 



(16) 



Using the eigenvalues of H given in equation(16) to equation(9), the term A^^- in equation(6) can be calculated 
(Ajj- is defined in equation(8)). From equation (16), it is seen that NSIs are present in the eigenvalues of H 
but those are of the order of a or a^. The oscillation probability Pii'e — ^ ^'ju) calculated upto order a^. 
As the non-zero terms in in equation(6) calculated using equation(7) is already at order a^, only terms 
which are zeroth order in a is considered in calculating A^^-. Hence, the oscillation probability, P(fe f/j,), 
in the presence of NSI upto order can be written as: 



Ac 



+ 



8S23C23 



23|6|2sin2' ^""31^'^ 



A'{1 - A') 

8.S23C23 



—-Re iba*) sin 



AE 
AmliA'L 



A'{1 - A') 



jzlTn (ba*) sin 



AE 



AE 



45^3 



{I- A 



l\2 



lap sin^ 



ArnJ^L 

AE 

Arni,L 
sm I I sm 



Amii(l- 


A')L^ 


AE 




Amii(l- 


A')L\^ 


AE 




'Am5x(l - 


A')L\ 


AE J 



(17) 



For NSI terms X = Y = Z = in equation (17), the probability expression reduces to that for the standard 
model interaction of neutrinos with matter. 



III. BIMAGIC CONDITIONS ON NEUTRINO ENERGY 



If we want the probability, Pivg ^n) to be independent of the CP violating phase 5 and ^13 upto order 
then we have to use the condition: 

™^=41i-i2i = o (18) 

AE 

One may note here that this corresponds to two different conditions for two different hierarchies of neutrino 
masses. Considering Cee = in A' and keeping e^^ and eer less than a as required by the perturbation theory 
one can see that this condition on neutrino energy is the same one as discussed earlier without NSIs [7]. For 
a given length L of the baseline above condition constrains the neutrino energy E as 

E = Amlil {±An-K/L + 2\/2GFne (1 + Cee)) ■ (19) 
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As long as e^e is unknown it seems that this magic energy E cannot be known. However, what is important 
in our work is to know the possible range of this magic energy depending on the presently allowed range of 
Eee which Is less than 4. Using eq.(18) in eq.(17) we get 




E (GeV) E (GeV) 



FIG. 1: Plots for probability of oscillation i/e — > Vy, versus neutrino energy {E) varying 5 over the entire range of 

to 2-n and 6\3 in the entire allowed range to 12° . 



In figure 1 we have shown the probability of oscillation Ve — ?• versus energy after varying Oiz in the 
entire allowed range of 0° — 12° and also varying 5 in the full range of — 27r for L = 2540 Km for both normal 
(NH) and inverted hierarchies (IH). We have made the plots for two NSIs (ce^ and e^r) but considering one 
at a time. The value of the NSI considered is = 0.024 and e^r = 0.024. From the plots we can see that 
the probability becomes almost independent of 6 and Oi^ at the magic energies. As for example, for n = 1 
in eq.(19) these energies are {E = 1.9(3.3) GeV for NH(III) respectively. Magic energies for higher n values 
gets smaller as seen from eq.(19) and the figure. 

Next we discuss the bimagic conditions on neutrino energy. As discussed in [7] the sensitivity of the hier- 
archy is maximum if one of the hierarchies (say NH or IH) obey the condition in cq.(18) for which probability 
is independent of ^13 and 6 whereas for the other hierarchy the probability has maximum dependence with 
^13 and 5 which can be achieved by imposing the condition 

Sin^=%^ = ±l. (21) 

These two conditions mentioned in (18 ) and (21) can be rewritten in two different ways: (a) conditions for 
IH with 6 and ^13 independence and NH with maximum dependence to those - which can be written as: 

»«=„,for.H (22) 
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^^l^\^'-\^'\)Ii = i^m-l/2).iornn (23) 



AEhc 

where n, m are integers and n > 0. 

(b) conditions for NH with 5 and ^13 independence and IH with maximum dependence to those - which 
can be written as: 



\^ml,\{l^\A'\)L 
AEhc 



= n'lT for NH (24) 



^--l(l + m)^ = (^._l/2).fo,IH (25) 



AEhc 

where n' and m' are integers and n' ^0 and m' > 0. 

Solving the equations (22) and (23) one gets the length of the baseline (L) as 

i(Km ) = l!i^^^i±l/^ « 16260.5 X (-"-+1/2) 



V2GFne{l + eee) ' p{gm/ Cc)(l + Cee) 

which implies 

pL(Km gm/cc) w 16260.5 x (n - m + l/2)/(l + Cee) (26) 
and for the inverted hierarchy the energy Eih with 5 and ^13 independence as 

F (CcV) ^ |Amii|(GeV^)L(Km) 

Similarly, solving the equations (24) and (25) one gets the length of the baseline (L') as 

pi' (Km gm/cc) w 16260.5 x (m' - n' - l/2)/(l + Cee) (28) 
and for the normal hierarchy the energy Enh with 6 and (^13 independence as 

F (CrV)- ^ |Amii|(GcV^)L-(Km) 

£;;vH(GeV) - 2^;,c(GeV Km) {n' +m'- 1/2) ^^^^ 

Firstly, for Cee = one can get one possible solution for common baseline i.e., L = L' to be about 2540 Km 
if the choices are made as follows: n = 1, m = 1 and n' = 1 and m' = 2 . The neutrino energy Ejh ~ 3.3 
GeV and EffH ~ 1-9 Gev. However, one could get more common solutions for bimagic baseline with L = L' 
for length by considering suitable choices of m, n, m' and n' for which n — m = m' — n' — 1; n — m could be 
1 or 4 resulting m L = V > 6000 Km. As for example, considering p w 4 gm/cc with n — m = 1 the length 
is about 6100 Km. 

As the present upper bound of some of the NSIs could be quite large a [13] which is not considered in 
our perturbative approach it is natural to ask what happen to such magic energies in the same 2540 Km 
baseline in those cases. We have checked numerically that even for highest allowed values of NSIs as for 
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examplc,for = 0.33 the Enh{Eih) « 2.03(3.26) GcV; for Cee = 4.0 the Enh{Eih) ~ 1.18(3.18) GeV. 
However, for e^r ^ 0.5 it is difficult to get bimagic energies although we have presented numerical analysis 
of discovery limits of various oscillation parameters for that also. For = 0.5 the bimagic energies are the 
Enh{Eih) ~ 2.02(3.57) GeV. It is important to note that all these bimagic energies are within 1-5 GeV 
which is the full neutrino energy range in our analysis. Interestingly, one can check here that the perturbative 
results for bimagic baseline length and the energies hold good even for higher values of egg >> a. Even for 
such large €ee it is possible to obtain bimagic energies within 1-5 GeV for the same baseline length of 2540 
Km. As for example, for Cee = 4 considering n = 3, m = l, m'=4 and n' = 1 gives L = L' w 2540 Km 
and EiH and Ej^h obtained from perturbative; approach are very near to the numerical values for bimagic 
energies mentioned above. As the bimagic energies are within 1-5 GeV, by considering this as the neutrino 
energy range in our numerical analysis, it may be expected to get better discovery limits to hierarchy, 9i3 
and CP violation for various choices of NSIs. 

Now, we can write down the probabilities at the particular energies Ejh and Enh according to the 
conditions discussed above. At E Eih for condition satisfying eq. (22) the probability is given by: 

P{v, ^ v,){IH) « sin^ (^|(^n + m - 0) (30) 

and for the condition satisfying eq. (23) the probability is given by: 

P(.e ^ ^,){NH) « ^|6pcos^ f; + ^l) + iShvM' 



|^'|2' ' \2\ 2) J (1- |A'|) 

8S23C23 



\A'\{l-\A'\) 



Re {ha*) cos^ ( ^ ('^ + ™~^)) + i^^*) sin ( tt ( n -|- m — 



(31) 



,2 V 2) j 1 ^ ' \ \ 2 

hk, E Ej^H for condition satisfying eq. (24) the probability is given by 

P{v, ^ v,){NH) « ^|6psin^ (| [n' + m' - 0) (32) 

and for the condition satisfying eq. (25) the probability is given by 



P{u, ^ u^){IH) 

8S23C23 
^\A'\{l + \A'\) 



Re (6a*)cos^ 



4si3 



n + m 



(l + l^'l) 
\-Iin {ba*) sin ( tt 



■ TO 



(33) 



One can see that at Ejh = 3.3 GeV from eqs. (30) and (31) in the case of c,.r = fe/j = case i.e., 
X = Y = a Oi3 also vanishes then a = and there is no difference in Pjh and Pnh- Same thing happens 
at Enh = 1.9 GeV as seen in eqs. (32) and (33). This means that for egr = Ee/* = case it is not so likely 
to get the hierarchy discovery limit at ^13 = 0. However, on the contrary in presence of these NSIs there is 
difference in Pjh and Pnh even for ^13 = 0. So in presence of NSIs like ggr and one could get discovery 
limit of hierarchy even at (?i3 = 0. However, in case of e^e if ^13 = then there is no difference between Pih 
and Pnh and so it is not likely to get discovery limit at 613 = 0. These features are found in our numerical 
analysis. 
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In probing other NSIs like e^^, e^,- and e^T-, oscillation channel z/,. — > i^^ is not appropriate. This follows 
from the probability in equation (17). For those NSIs considering the disappearance channel (i/^ — >■ f^) is 
appropriate. One cannot get any condition on neutrino energy in general to remove the dependence on S in 
the oscillation probability for this channel. However, it is found that upto the order a without imposing any 
condition on neutrino energy this oscillation probability is already independent of d as shown below [15] . 

, , .99.9 /^Tfi\-t L .9 9 9 ^TTloT L . ^TTlo-t L 

„ 9 9 r/- 2 2 \ / \ . r. / ^^ AmoiA'L . AmoiL 
+ 2C23S23 [(C23 - S23) (e^ - err) - 4c23S23i?e(e^i-)J — sm 

- 8C23S23 (c23 - 5^3) [C23S23 (e^p - err) + (c^3 " sIs) ^e(e^^)] A' sin^ (34) 

To get the sensitivity of NSIs like e^^, e^^ and c^t, this disappearance channel is appropriate but it is not 
much suitable for finding sensitivity to 5 or discovery limits for CP violation. As the probability upto order 
a is already independent of S here one does not require the bimagic conditions and as such there is no 
restriction on the length of the baseline and neutrino energy. We have not done this analysis separately but 
sensitivity of some of the above-mentioned NSI parameters have been discussed in the disappearance channel 
(i^M f/x) in [16]. 



IV. NUMERICAL SIMULATION 



As an outcome of the bimagic energy conditions in presence of NSIs the length of the baseline is ~ 2540 
Km. To study the oscillation of Ug — )• i/^ we have considered the experimental set-up and the detector 
characteristics as discussed in [7] for a running time of 2.5 year. We consider the neutrino factory having 
5 X 10^^ muon decays per year with parent muon of energy 5 GeV and the magnetized totally active scintillator 
detector of 25 kt mass with energy threshold of 1 GeV. The numerical simulation has been done by using 
GLoBES [17]. In presenting the discovery limits of hierarchy, sin^ ^13 and CP violation in bimagic baseline, 
highest possible values for the NSIs have been considered. However, when we observe that the discovery 
limits are either covering the entire allowed region or not at all obtainable then we refrain from presenting 
those figures and instead we present discovery limits for some lower values of NSIs. 




le-05 0.0001 ^ 0.001 0.01 1S-05 0.0001 ^ 0.001 0.01 le^K 0.0001 ^ 0.001 0.01 

sin (trae) sin 613 (trae) sin 613 (trae) 



FIG. 2: The 3a contours showing discovery limit of hierarchy for different NSIs Cer, ee^i and e, 



In figure 2 we have shown for what value of 6 and ^la at 3 a confidence level one can identify the specific 
hierarchy which could be either normal or inverted. For Cer = 3 and = 0.33 (one at a time) NH could 
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be discovered at any vahie of 6is irrespective of any specific value of S whereas for IH nowhere it is found to 
be discovered. So in our figures we have chosen some lower values of these two NSIs. Considering e,,f^ =0.1 
and = 0.1 (one at a time) we have shown in figure 2 the discovery limit of hierarchy. Prom these figures 
for favorable values of S one could identify the inverted hierarchy of nature for Cer at sin^ 613 as small as 
4 X 10~^ and for Cef^ at sin^ ^13 as small as about 2.5 x 10~^ . For normal hierarchy however, it is found from 
the figures that for eer =0.1, only for 6 in the range of 37r/4 to 57r/4 and sin^ di^ 5 x lO^"' one could reach 
the discovery limit. For other values of 6 normal hierarchy can be identified for any value of 613 including 
the zero value. Similarly, for ee^u = 0.1, only for 6 in the range of about 7r/4 to 57r/4 and sin^ ^13 ^ 2 x 10~^ 
one could reach the discovery limit. Here also for other values of i5 normal hierarchy can be identified for any 
values of 6*13 including the zero value. For e^e = 4 the normal hierarchy can be identified at sin^ ^13 ^ 10"'' 
and the inverted hierachy can be identified at sin^ ^13 ^ 10~^. 




sin^ (true) sin^ (true) sin^ 613 (true) 

FIG. 3: The 3a contours showing discovery hmits of ^13 for different NSI's and te 




FIG. 4: The 3a discovery limits for the CP violating phase S for different NSI's eer, ee^ and eee- 



We have shown the discovery limits of ^13 in figure 3 at 3a confidence level. From figure 3 it can be seen 
that for NSI 

^en 0.33 the discovery limits for sin'^ ^13 could be as low as sin'^ O13 ^ 5.0 x 10 ^ and for 

Cer 3.0 the limit could be as low as 2.0 x 10"'^ for both the hierarchies. However, for e,;e ^ 4 this limit 
improves for NH and can be as low as sin^ ^13 ^ 2 x 10"''', but for IH it could be as low as sin^ 613 ^ 1.5 x 10^"^. 
In figure 4 we have shown the discovery limit of CP violation for different NSI at 3a confidence level. The 
discovery limit for CP violating region is possible for Cer = 3.0 at sin^ ^13 > 2.5 x 10~^ for 37r/8 < ^ < 37r/4 
for NH. But for IH it is very difficult to obtain any discovery limit. However, for lower values of egr discovery 
limits could be easily obtained. For e^^ = 0.33 the discovery limit of CP violating region for NH is found for 
sin^ 6'i3 ^ 2 X 10-3 with Tr/8^6<, Ttt/S and also for sin^ 6*13 ^ 8 x lO"'' with 97r/8 ^ 6 <, Ibn/S. In the case 
of IH for the same value of the discovery limits are found for sin^ ^13 ^ 4 x 10"^ with 7r/8 ^ S ^ 37r/4 
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and also for sin^ ^13 ^ 4 x 10 with 57r/4 ^ S ^ 157r/8. For higher value of e^e = 4 it is not possible to get 
discovery limit for CP violating region. For e^,, = 0.4 the discovery limit of CP violating region for NH is 
found for sin^ 6*13 ^ 7 x 10"^ with n/A <, 6 <, tt/2 and also for sin^ 6I13 ^ 6 x 10"^ with 5tt/4: <, 6 <> 7tt/A. 
In the case of IH for the same value of Cee the discovery limit are found for sin^ ^13 ^ 9 x 10~^ with 
7r/8 ^ 5 ^ 7r/2 and also for sin^ 613 k, 2 x IQ-^ with 57r/4 ^5^ 77r/4. 




FIG. 5: Discovery limits of NSI (ser) for different fixed values of ^13 and S considering muon energy 5 GeV. 



Although wc find good discovery limits for hierarchy, ^13 and CP violation in bimagic baseline even 
in presence of NSI, however, to get good discovery limits of NSI, the neutrino energy around 5 GeV (as 

required by magic energy conditions) is not appropriate. One can see from the expression of P{i>e — > I'p,) 
in Eq. (17) that the NSI terms are energy independent whereas the terms containing only vacuum mixing 
parameters are suppressed by neutrino energy. This feature is present irrespective of specific channel for 
neutrino oscillation. Naturally for higher energy the relative effect of NSI parameters are enhanced with 
respect to vacuum neutrino mixing parameter and one might expect to get better discovery liimits of NSIs. 
Now considering 5 GeV as the maximiim neutrino energy, we have presented the discovery limits of some 
NSIs - €eT, iefj. and Cee foi" various fixed values of sin^ 613 and 6 in figures 5, 6 and 7 respectively . We can see 
from these figures that the discovery limits of Cgr and are as low as (w 0.015) for either of the heirarchy 
at 3(7 confidence level. For Cee the limit is as low as at the order of (~ 10~^) . But at higher neutrino energy 
say for 50 GeV from figures 8 and 9 one can see that for IH the discovery limit of egr can be as low as 
3 X 10~^ and that of eg/* could be as low as 7 x 10"'*. Similarly, for NH the discovery limit of eer is as low as 
0.01 and for Cgfj, is as low as 0.002. For the case of Cee, from figure 10 we can see that the discovery limit of 
the NSI (cee) is not so good and could be as low as of the order of 10^^. However, the overall probability of 
oscillation is suppressed with the increase in neutrino energy. Naturally it is expected that just increasing 
energy one may not keep getting better NSI discovery limits. In fact, we have checked at neutrino energy 
above 60 GeV there is insignificant improvement in discovery limits of and Cer in 2540 Km baseline. 
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FIG. 6: Discovery limits of NSI (se/i) for different fixed values of 613 and 6 considering muon energy to be 5 GeV. 




PIG. 7: Discovery limits of NSI (eee) for different fixed values of ^13 and S considering muon energy 5 GeV. 




FIG. 9: Discovery limits of NSI (se/j) for different fixed values of ^13 and 5 considering muon energy 50 Gev. 
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FIG. 10: Discovery limits of NSI (tee) for different fixed values of ^13 and S considering muon energy 50 GeV. 

V. CONCLUSION 

It is found that for getting good discovery limits for hierarcliy, sin^ ^13 and CP violation particularly in 
the Uf. — )• oscillation channel, 2540 Km baseline is suitable even when NSI of neutrinos with matter are 
present. This is because the bimagic energies Ejh and Enh lie within specific energy range, which is 1-5 
GeV for this baseline even in presence of NSIs with their lower or higher allowed values ( except for e^r ^ 0.5) 
and this neutrino energy range has been chosen in our analysis with NSIs. It is important to note that this 
energy range is also suitable for no NSIs as in that case also bimagic energies are within 1-5 GeV [7]. 

To show what could be the utmost effect to the discovery limits corresponding to no-NSI case, we have 
considered highest possible values as obtained in the model independent case [13]. However, for model 
dependent cases [8-10] these bounds are in general, more stringent. 

The discovery limits of hierarchy actually improves in presence of NSIs. Even one could get discovery 
limits at 61^ = for e,,^ and e^r which in absence of those NSIs are not expected. This is due to the fact 
that at bimagic energies the Pjh and Pnh are unequal even at 813 = in presence of those NSIs. This 
does not occur for Cee- In this case, as for example, for Cee = 4 the hierarchy discovery limits could be 
obtained at as low as sin^ 613 ^ 10"''. Considering highest possible allowed values of Cer, ^ep. and e^e we 
find that the discovery limits of sin^ ^13 could be as low as 2 x 10""^, 6 x lO^"' and 2 x 10~^ respectively 
for normal hierarchy and as low as 2.8 x 10"^, 7 x 10""* and 1.5 x 10""^ respectively for inverted hierarchy. 
Considering favorable values of S the discovery limits of CP violation are possible at following sin^ ^13 values. 
For = 3 the discovery limits of CP violation could be possible for high value of sin^ O13 at about 0.025 
for normal hierarchy only. For inverted hierarchy it is not possible. For = 0.33 the discovery limits 
of CP violation could be obtained for sin^ ^13 as low as 10~^ for normal hierarchy and at about 4 x 10~^ 
for inverted hierarchy. For Cee = 4.0 discovery limits of CP violation cannot be obtained. However, for 
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lower values of both eer and one could get discovery limits of CP violation at some sin^ ^13 values. The 
discovery limits of NSIs could be improved if we consider neutrino energy upto 50 GeV and it could be as 
small as 10~^ for e^ij, and and could be as small as 10~^ for Cee- These NSI discovery limits essentially 
would give the upper bound on the respective parameters if they are not discovered. 

It is interesting to note that there are other bimagic baselines with length greater than 6000 Km apart from 
2540 Km as discussed before. One may explore the discovery limits of various vacuum neutrino oscillation 
parameters using those baselines also. 
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